An Upper Limit on the Granularity of the Local Stellar Halo 
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I use the statistical properties of 4588 nearby halo stars to show that if the 



\ local stellar halo is composed of cold streams from disrupted dwarf satellites as 



is predicted by popular cosmologies, then at 95% confidence there are at least 



\ 1350 such streams. Moreover, no stream can contain more than 2.3% of the local 

^ ' stars, while average streams must contain less than 0.08%. 
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1. Introduction 

In the conflict between the two competing models for the formation of the Milky Way, 
the Searle & Zinn (1978) scenario of cannibalistic assimilation of pre-existing protogalaxies 
has generally been gaining ground on the earlier monolithic collapse model of Eggen, Lynden- 
Bell & Sandage (1962). In part this is because the former fits naturally into the picture of 
hierarchical structure formation predicted by cold dark matter models, which have so many 
successes on larger scales. Even more compelling is the direct observational evidence. The 
best example is the Sgr dwarf galaxy (Ibata, Gilmore & Irwin 1994), whose disruption is 
taking place right before our eyes. However, while it is especially striking (Majewski et al. 
2003 and references therein), the Sgr dwarf is by no means unique. Palomar 5 (Odenkirchen 
et al. 2001; Rockosi et al. 2002) and possibly the Ursa Minor dwarf spheroidal (Palma 2003) 
show extended tidal tails as evidence of ongoing disruption. 



However, these disintegrating satellites are all in the outer halo. While individual moving 
groups have been reported in the local stellar halo for several decades (Eggen 1965), it is not 
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clear whether these associations are an exotic froth on an underlying well-mixed homogeneous 
stellar population, or whether the local stellar halo is basically granular, being composed of 
a number of streams, each tightly clumped in velocity space but so extended in physical 
space that they cannot be distinguished photometrically from the myriad of other streams. 

Numerical simulations demonstrate that kinematically cold, initially spherical systems 
will form cold streams when tidally disrupted (McGlynn 1990), and such disruption has been 
studied in considerable detail (Johnston 1998). Bullock, Kravtsov & Weinberg (2001) sim- 
ulated the formation of a network of streams within the framework of hierarchical structure 
formation and showed that halo tracers such as RR Lyrae stars would exhibit dramatic vari- 
ations in surface density projected on the sky. While these variations are most pronounced 
at large distances, they are also clearly discernible in the simulated solar neighborhood. Sim- 
ulations of satellite disruption by Helmi & White (1999) also support the granular picture. 
They specifically predict a local halo composed of 300-500 streams, each with a velocity 
dispersion of order 5kms~^. 

Here I show that if the local stellar halo is composed of such streams, they are extremely 
numerous, with typical streams each contributing less than 0.08% to the local stellar halo. 

2. Limits on Granularity 

In Gould (2003), I analyzed a sample of 4588 halo stars selected from the revised NLTT 
catalog (rNLTT) (Sahm & Gould 2003; Gould & Sahm 2003) by means of an optical-infrared 
reduced proper motion (RPM) diagram (Salim & Gould 2002). In particular, I found that 
five of the nine parameters describing the local halo velocity ellipsoid were consistent with 
zero. These were the mean radial and vertical motions of the halo relative to the local 
standard of rest (LSR), Vi — 1.4±2.2kms~^ and V3 = 1.8 it 2.4kms~^, and the normalized 
off-diagonal components components of the velocity dispersion tensor, C12 = 0.024 ± 0.014, 
Ci3 = 0.005 ± 0.023, and C23 = —0.004 ± 0.026. (Of the remaining four velocity-ellipsoid 
parameters, the fit values of the three velocity dispersions were found to be compatible with 
previous determinations and the motion in the tangential direction relative to the Sun was 
fixed as part of the analysis at the previously determined value of U2 = — 216.6kms~^.) 

There are many possible effects that could have driven these five numbers away from 
zero, which is the value predicted by the simplest axisymmetric Galactic model. For example, 
if the LSR were moving inward or outward relative to the Galactic potential, the reflex motion 
of the halo would have driven Vi away from zero. Or if the halo had been rotating around the 
2-axis, this would have shown up in V3. Or if the dark halo were significantly triaxial with 
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principle axes misaligned relative to standard Galactic coordinates, this would have driven 
Ci2 to significant non-zero values. However, and this is the key point, there are no systematic 
effects that could have driven these parameters systematically toward zero if that were not 
their actual value. In particular, the likelihood routine actually fit for motions relative to the 
Sun, which were found to differ from zero by 5.2 and 2.3 sigma. It was only after the known 
motion of the LSR was removed that they became consistent with zero. Since these five 
parameters agree with the predictions of the simplest model, their statistical fluctuations, as 
codified by x^, can be used to test the statistical properties of the underlying sample. Note 
that Xnom — 3-97 for 5 degrees of freedom (dof). 

Let us now suppose that the local halo were composed of 459 streams, each with 10 stars 
and each with an internal velocity dispersion well below the 100 km s"^ to 170kms~^ values 
measured for the three components of the halo. The 10 stars would then have had essentially 
identical kinematics, and I would therefore have overcounted the number of independent 
measurements by a factor of 10. The best fit parameters for the halo that I derived would 
have been correct, but the errors would have been underestimated by a factor of 10^/^. Hence, 
would have been x^ruc = Xnom/10 = 0.397. This scenario can be immediately ruled out 
because the probability of finding such a low ior 5 dof is less than 0.5%. 

More generally, one can imagine that the 4588 stars in the sample are divided into 
streams, with the ith stream having rrii members, > 1. Then would have been 
overestimated by ^^771^/4588. Hence, 

2 4588 2 

Xtrue 2 Xnom 

I now demand that Xtrue > 1-145, so that it does not fall within the 5% lowest part of the 
X^ distribution for 5 dof. I thereby obtain a limit, 

J2 m\ = ^ J] mi < 15908, J] = 4588. (2) 

i=i '^'^t^e i=i i=i 

If the nii were real variables, the minimum value of N that would saturate this constraint 
would be = 4588715908 = 1323, with = 15908/4588 = 3.47. Since the must be 
integers, the minimum is N — 1350, in which case there would be 816 streams with 3 stars 
and 535 with 4 stars. That is, each stream would have about 0.08% of the total local halo 
population. The maximum number of stars that can be in any stream is 106. In this case, 
all the other 4483 "streams" would have only one star each. 
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3. Discussion 

If the statistical errors in Gould (2003) had been overestimated, this would have arti- 
ficially decreased ^"^^ so made the limits derived in the present Letter appear stronger 
than they actually are. Is such an error overestimation a possibility? No. The Gould (2003) 
errors are strongly dominated by Poisson statistics, not measurement errors, so that even if 
the measurement errors were set to zero in the analysis, the errors on the derived parameters 
would not be reduced by more than a few percent (Popowski & Gould 1998). 

Suppose that the Gould (2003) halo sample were contaminated at a fractional level / 
by another population (e.g. thick disk stars) whose five velocity-ellipsoid parameters were 
also zero. Suppose that the halo sample were clumpy but the other population was not. The 
lower hmit derived in § 2 should then be decreased by a fraction /(r — l)/(r — /) ~ 0.7/, 
where r = Xnom/x^(5%) = 3.47. Prom figure 3 of Salim & Gould (2003), I estimate / ~ 2%, 
so this effect is negligibly small. 

Gan the technique presented here be extended to probe more effectively for granularity 
in the stellar halo? The rNLTT covers only 44% of the sky and moreover is derived from a 
proper-motion selected catalog /i > 180masyr~^. If the catalog were extended to most of 
the rest of the sky (see Salim & Gould 2003), the number of stars could be doubled. By 
comparing first epoch Palomar Observatory Sky Survey (POSS I) plates with Early Release 
Data from the Sloan Digital Sky Survey (SDSS), Digby et al. (2003) were able to identify 
halo stars from a RPM diagram for proper motion stars down to 40masyr~^. This extends 
the distance out to which they can be found by a factor 4.5 relative to rNLTT and therefore 
increases the volume by a factor 90. Because halo stars are distributed roughly uniformly in 
space, the sample size will increase in proportion to the volume. Since SDSS will eventually 
cover 1/4 of the sky, the prospects for increasing the sample by a factor ~ 50 are good. 
Future surveys may do even better. If the five velocity ellipsoid parameters continue to 
be consistent with zero, it will be possible to push down the limits on halo granularity by 
roughly the same factor. 

But what if some of these parameters begin to deviate significantly from zero? How can 
one tell, within the framework of this approach, whether these formally significant detections 
refiect the real properties of the halo velocity ellipsoid or an underestimation of the velocity- 
ellipsoid errors due to dumpiness of the stellar halo? 

Two of the five parameters, Vi and Ci2 are much more likely to really be non-zero than 
the other three. If the halo is triaxial and one of its axes is not perfectly aligned with the Sun, 
then this should induce both inward (or outward) net motion of the LSR (Vi) and a vertex 
deviation in the halo (£12). If one or both of these parameters begins to differ significantly 
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from zero, but the remaining three remain consistent with zero, then the same test developed 
here can be used to constrain granularity, except that it would have to be restricted to the 
three remaining parameters, V3, C13, and C23. On the other hand, if halo granularity is causing 
the errors to be underestimated, it should affect both groups of parameters equally. 

While the statistical power of a test for 3 dof is somewhat weaker than for 5 dof, this 
would be more than compensated by the increase in star counts. Thus, there are reasonable 
prospects for simultaneously carrying out high-precision tests of Galactic structure using two 
of the five parameters while tightly constraining or precisely measuring the degree of halo 
granularity using the other three. 

The limits derived here are in strong formal conflict with the predictions of Helmi & 
White (1999). No formal comparison can be made with the predictions of Bullock, Kravtsov 
& Weinberg (2001) because they do not report on the quantity measured here, the total 
number of cold streams in the local halo. It is nevertheless possible that a prediction for this 
quantity could be extracted from their simulation. 

However, this formal conflict does not necessarily contradict the picture of hierarchical 
structure formation on galaxy scales. The simulations conducted to date have considered 
one or a few realizations without an attempt to tune parameters to fit any quantitative 
measurement of streams. With the advent of the limits presented here, as well as the 
prospect of much better limits and/or measurements in the reasonably near future, it should 
be possible to quantitatively constrain such scenarios or, perhaps, to rule them out. 

I thank Scott Tremaine and Scott Gaudi for stimulating discussions. This work was 
supported by grant AST 02-01266 from the NSF and by JPL contract 1226901. 
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